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v^ Abstract 
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We provide analytical results for the probability distribution of a family of wavefunc- 
_^ tions of a quantum mechanics model of commuting matrices in the large-TV limit. 

These wavefunctions describe the strong coupling limit of 1/8 BPS states of M = 4 
Oh 

supersymmetric Yang-Mills theory. In the large- N limit, they should be dual to clas- 
sical solutions of type IIB supergravity that asymptotically approach AdSs x S 5 . Each 
probability distribution can be described as the partition function of a matrix model 
(different wavefunctions correspond to different matrix model potentials) which we 

study by means of a saddle point approximation. These saddle point solutions are 
OO 
P^j given in terms of (five-dimensional) hypersurfaces supporting density distributions of 

^O eigenvalues. 
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1. Introduction 

Since the discovery of the anti-de Sitter space/conformal field theory (AdS/CFT) correspondence 
|1|, significant progress has been made in our understanding of strongly coupled gauge theory 
phenomena. So far, the example that has been studied most extensively is the correspondence 
between SU(A r ) maximal AA = 4 supersymmetric Yang-Mills (SYM) theory in four dimensions and 
type IIB superstring theory on the ten-dimensional AdSs x S 5 background. The correspondence 
in its strongest form claims full dynamical agreement between both theories at the quantum level. 
This is certainly hard to verify and one therefore seeks limits in which it is possible to perform 
tests explicitly. One very interesting limit is the planar or large- A limit of the gauge theory. In 
this limit, N = 4 SYM theory is believed to be equivalent to free (i.e. genus-zero) string theory on 
AdSs x S 5 . If one in addition assumes strong 't Hooft coupling (or equivalently, large curvature 
radius of AdS), then the gauge theory is believed to be described by classical supergravity. 

One may adopt the point of view of taking the AdS/CFT correspondence as an approach 
towards defining a theory of quantum gravity. This then naturally leads to the question of how 
geometrical information emerges in the strong coupling limit of the gauge theory. Notice that 
gravity is not apparent in the Lagrangian description of quantum field theory and in this sense 
it can be thought of as 'emergent'. Moreover, the gravity side of the correspondence is higher 
dimensional, therefore one should also be able to understand how excitations are localised in the 
dual extra dimensions by field theoretical means. 

To shed light on the problem of emergent geometry, one needs to study strongly coupled 
gauge theory, which on general grounds, is extremely hard. To tackle this problem, Berenstein |2| 
proposed a truncation of N = 4 SYM to a quantum mechanical problem of commuting matrices 
by compactifying the theory on a three-sphere (i.e. he considered J\f = 4 SYM theory on R, x 5 3 ). 
This compactification provides a natural infra-red regulator. Moreover, upon expanding all the 
fields of N = 4 SYM theory in spherical harmonics on the three-sphere, the action is truncated 
to obtain a quantum mechanics Hamiltonian of six Hermitian matrices. This truncation of the 
degrees of freedom to commuting matrices is a good approximation in the strong coupling limit. 1 

The notion of emergent geometry in the strong coupling limit of more general matrix models was studied in 3 . 



The eigenstates of the model with six commuting matrices are conjectured to describe 1/8 BPS 
states of M = 4 SYM theory at strong coupling. In the large- N limit, these are dual to classical 
solutions of type IIB supergravity that asymptotically approach AdSs x S* 5 . Similar models have 
been proposed to describe BPS states in orbifolds of J\f = 4 SYM theory |4| and certain J\f = 1 
superconformal field theories |5|. 

A Gaussian wavefunction of the eigenvalues of the commuting matrices was shown to be the 
exact ground state of the quantum mechanical Hamiltonian 2 1 . The product of this ground state 
by a holomorphic function of the matrices' eigenvalues is an eigenfunction of the Hamiltonian, 
as well, at least to a good approximation in the large- N limit |5,6|. Now, having this family of 
quantum mechanical eigenstates, it would be desirable to characterise their typical or most likely 
distribution of eigenvalues. The probability distribution of each of the wavefunctions can be seen 
as the partition function of a given matrix model. For the ground state, the resulting partition 
function would be that of a model with a quadratic potential and a generalised Vandermonde 
repulsion. For the 'excited' wavefunctions, other terms are added to the potential of the model. 
In the large- N limit, all these partition functions will be dominated by their saddle points. Al- 
though only for the ground state the saddle point equations have been solved exactly so far |2,7|, 
the result is very compelling: The saddle point configuration is a uniform distribution of eigen- 
values supported on a five-sphere embedded into R 6 . Just to remind the reader, the ground state 
wavefunction must be identified with the dual of the AdSs x S 5 background. Subsequent studies 
revealed that it is possible to confer an explicit geometrical interpretation to the five-sphere of 
eigenvalues |7|. In order to see if one can push this identification further, it would be necessary 
to compute the partition functions associated with other eigenstates of the quantum mechanical 
Hamiltonian. The simplest cases one could start from to look at are 1/2 BPS eigenstates. For 
these states the AdS/CFT dictionary has been studied thoroughly (see e.g. |8|) and it is known 
how to relate them to a family of 1/2 BPS supergravity solutions that asymptotically approach 
AdSs x S , found by Lin, Lunin &; Maldacena (LLM) 1 9 2 , Unfortunately, solving the corres- 
ponding partition functions in the quantum mechanics model with six commuting matrices is a 
very difficult problem, even in the saddle point approximation. The only reported results in this 
direction are in a series of articles that study some of these wavefunctions numerically for finite 
N [6,12,13]. 

In this paper, we will extend the analytical result that is known for the ground state wave- 
function to excited wavefunctions. Specifically, we will develop a perturbative method that allows 
for an analytical treatment of the saddle point equations of a family of wavefunctions. We will 
present analytical results for monomial and logarithmic potentials in the corresponding matrix 
model Hamiltonian. Wavefunctions for degree pgN monomial potentials are expected to be 
in correspondence with LLM geometries obtained from a simply connected droplet possessing a 
non- vanishing p th harmonic moment |10|. Wavefunctions for logarithmic potentials should cor- 
respond to annular LLM geometries where the inner radius is determined by the strength of 
the logarithmic potential |2|. Our results should form a starting point for extracting geometric 
information directly from the gauge theory. We will comment on this issue in the conclusions. 

This paper is organised a follows. We will first give a brief review of the background material. 
In Section 3., we then provide our perturbative approach and present analytical solutions for 
monomial and logarithmic potentials. We shall also compare our results against the numerics. In 
Section 4. , we will conclude and give an outlook of open problems. Finally, several appendices 
collect useful definitions and details of our derivations. 

2 An alternative approach to recast geometrical information of LLM solutions was used in 10, 11 . 



2. Quantum mechanics of commuting matrices 

The system we are going to be dealing with is a particular matrix quantum mechanics model of 
six commuting Hermitian matrices. Let X = (X , . . . , X 6 ) be the six Hermitian matrices, they 
are therefore subject to the constraint 

X AX = . (2.1) 

The Hamiltonian we are interested in is 

h cX = |tr(n-n + x-x) , (2.2) 

where II is conjugate to X and 'cl' refers to classical. The system has a gauge invariance, where 
one acts by conjugation: X \-t g~ 1 Xg for g 6 SU(iV). Because of the constraint (2.1), one can use 
this SU(iV) action to diagonalise all six matrices X simultaneously. Let us denote the eigenvalues 
of X by Xi = (xj, . . . , x\ ) E R 6 for i, j, . . . = 1, . . . , N. Having diagonalised all the matrices, we 
have fixed a gauge. However, there are still residual gauge transformations which permute the 
eigenvalues X{. Therefore, the corresponding wavefunctions will eventually be symmetric under 
the exchange Xi o Xj for all i and j. 

As shown by Berenstein 2 1 , this system can be obtained as a truncation of A/ - = 4 SYM theory 
onEx5 3 , to the s-wave modes of its six scalar fields, to describe gauge invariant 1/8 BPS states 
in the strong coupling limit. 

Having reduced the dynamics of these six N xN matrices to the dynamics of their eigenvalues, 
the system can thus be interpreted as set of N bosons on a space with six dimensions. If we treat 
the system classically, we can use a diagonal ansatz to find solutions of the dynamical system. 
Under these assumptions, we find N free harmonic oscillators in six dimensions, which should be 
treated as N identical particles (bosons) on a six-dimensional harmonic oscillator. 

Quantum mechanically, we cannot do that immediately. This is due to a certain measure 
factor that arises from the volume of the gauge orbit, and which affects the dynamics of the 
system. This measure factor was computed in |2| and it is given by 

I? := II l^- f jf • ( 2 - 3 ) 

l<i<j<N 

The resulting quantum Hamiltonian is therefore 

- / 1 1 \ 

H = E^-^^VV^-I^J. (2-4) 

2.1. Wavefunctions 

The main object of study in this paper will be some wavefunctions of the Hamiltonian (2.4). The 
presence of the Vandermonde measure factor makes the corresponding Schrodinger problem very 
difficult. Notice nonetheless that the rather simple wavefunction 




ip (xi,...,x N ) = expl --)\xi\ 2 ) (2.5) 



is an exact wavefunction of H. In fact, it is the ground state wavefunction, 

N 
H^ = E ijo , with E = —(N-1) + 3N. (2.6) 

3 



The measure factor fj,, which will appear in the probability density distribution, can be ab- 
sorbed into the wavefunction ip by a similarity transformation 

if, i-> ip := pip . (2.7) 

In the following, we shall concern ourselves with the re-scaled wavefunctions ip only. 

If we square ip, we get a probability density distribution on the phase space of the N particles. 
For the ground state ipo, this is given by 

\ipo\ 2 = /i 2 exp -^|fi| 2 = exp -^|xj| 2 + ^ \og\xi~Xj\ 2 . (2.8) 

V i=l / \ 8=1 l<i<j<N I 



If we set J??o := — log \ipo\ 2 , then 

Jo 



f j jjd 6 Xi] exp(-(3J4? ) , with = 1 (2.9) 

can be interpreted as the partition function of a gas of particles in a confining external quadratic 
potential, Yli \%i\ 2 , together with a logarithmic repulsion term, — J2i<j 1°& Wi ~ Xj\ 2 i between the 
particles in six dimensions. 3 

Less is known about the exact excited wavefunctions for the Hamiltonian H. However, in |2,6| 
it was shown that for 1/2 BPS states, the wavefunction 

-0/(fi,...,£jv) := exp{^2f(z i )\ip (x 1 ,...,x N ) (2.10) 

is an approximate eigenfunction of H in the thermodynamic limit N — >• oo, provided / = f(zi) is 
a holomorphic function of z% ■= xf -+- ixf. The partition function in this case is then given by 

3T f = f (fld 6 Xi j exp(-/3Jt? f ) , with = 1, (2.11a) 

where 

N N 

jr f := -logl^/l 2 = ^|f i | 2 -2^1He/(z i )- ^ log |f; - x,f . (2.11b) 

i=l i=l l<i<j<N 

Here, '*He' denotes the real part. 

In this work, we shall be interested in the large- iV limit of the 'matrix model' partition function 
(2.11a). In this limit, the bosons will form some type of distribution density p on the phase space 
of a single particle (density of eigenvalues) . The goal for us is then to determine the shape of the 
density p using a saddle point approximation. 

2.2. Large-iV limit 



Next we wish to perform the thermodynamic limit of the Hamiltonian (2.11b). When taking the 
limit iV — t- oo, we may trade the sums in Jfff for integrals at the expense of introducing a density 
p = p{x) (which is constrained to be non- negative) , 



J2 -> fd 6 xp(x) and ^ ->• ^ f d 6 xd 6 x'p(x)p(x') , (2.12a) 

*=1 l<i<j<N 

This logarithmic repulsion generalises the Vandermonde repulsion of eigenvalues in matrix models of 14 . 



and which is subject to the normalisation 



d 6 xp(x) = N . 



Therefore, (2.11b) becomes 

jf f = [d 6 xp(x)\x\ 2 -2 J d*xp(x)V\ef(z) 



d x d x p(x)p(x ) log \x — x 



•/i 2 



(2.12b) 



(2.13) 



Notice that the constraint (2.12b) might be added to (2.13) by using a Lagrange multiplier A. 

In the large-iV limit, the partition function (2.11a) will be dominated by its saddle point. 
Then, the most likely density p can then be obtained by extremising the Hamiltonian Jiff. Spe- 
cifically, the variation of Jtf with respect to p yields 



A + \x\ 2 



2«Ke/(z)- d b xp{x)log\x-x 



J\2 



o. 



(2.14) 



where the constraint (2.12b) is enforced by A. This is the integral equation that determines the 
density p. Upon acting with A 3 on this equation, where A is the Laplacian on R 6 , one quickly 
realises that p cannot be an ordinary function but must be of distributional support |2|, 4 Thus, 
the integral equation (2.14) can only hold in a suitable region of R 6 . In particular, for the ground 
state where / = 0, the density p is uniformly supported on a five-sphere S 5 C R 6 |2,7|. 



S 3 : (a, £i, 6) 




Figure 1: Five-sphere described by the spherical coordinates (2.15) (for fixed r): The hemisphere 
is parametrised by ui = (0, <f>) while the three-spheres are parametrised by (a, ^1,^2)- Notice that 
the radii of the three-spheres shrink to zero as 9 approaches zero. 

Instead of Cartesian coordinates (x , . . . , x 6 ), we will find it fruitful to make use of the following 
spherical parametrisation of R 6 \ {0}: 



r sin 6 sin a sin £1 



r sin 9 sin a cos £1 



r sin 9 cos a cos £2 > x = r cos 9 cos 4> > 



x 
x 6 



r sin 9 cos a sin £2 , 



(2.15) 



r cos v sin 1 



with 9, a G [0, 7r/2] and (ft, £12 £ [0, 27r). Notice that with this choice of coordinates, we have 

(2.16) 



x 5 + ix 6 



r cos u e 



To proceed we should make an ansatz for the density p. We will take it to be supported on a 
hypersurface in R 6 . Then, since in all the cases we consider the potential deformation / = f(z) 



4 In this argument, one uses the fact that AiJ log \x — 3f\ ~ <r (x — x"). 



depends only on the holomorphic coordinates z, we propose the following ansatz for the density 
P- 

P& = 2^,0)5 6{r ~ m ^ ^ *> ■ (2 - 17) 

Here, p = p(9, 4>) and r = f (0, (/>) are non-negative functions which depend only on and <j) and 
5 indicates Dirac's delta function. The problem of finding the density p has thus been translated 
into finding the functions p and f. Notice that the constraint (2.12b) then becomes 

/ i 7r/2 f2TT 

/ d0 / dcp cos9 sin 3 0,3(0,0) = 1. (2.18) 

Jo Jo 

In what follows, we shall simplify notation and make use the abbreviations ui := (9,4>) to- 
gether with 

/ 1 7r/2 /*27T 

dfi 2 := / d0 / d0cos0sin 3 0, f f f 

Jo Jo =► / dQ 5 := / dQ 2 / dfi 3 • (2.19) 

/ 1 7r/2 / 1 27T /*27T " J J 

d^3 := / da d£i / d^2 sin a cos a , 
Jo Jo Jo 



Next we wish to substitute the ansatz (2.17) into (2.13) and derive the corresponding equations 
of motion. To this end, we need the expression 



I -> ^|2 
\X — X \ 



|x1 = , H = r 2 (u;) + r 2 (u')-2r(u;)f(u;')cos<p, (2.20) 

|z'|=r(u/) 



where (p is the angle ^(]§r, t§tt) between r|r and t^t on the five-sphere defined by 



cose/? := cos cos 0' cos (^ — </>') 

+ sin sin 0' [cos a cos a' cos(£i — £i) + sin a sin a' cos(^2 — £2)] ■ 

Therefore, we obtain for (2.13) (including the Lagrange multiplier term) 
Jf f = A( /d0 2 p(w)-l) +N fdn 2 p(u) [r 2 (u)-2 9\ef(z)] 

- ^j /dft 5 /W 5 pH/V) log [f 2 M + r V) - 2 f (w)f (a/) cos p] , 



(2.21) 



(2.22) 



where z = z(oj) = f(io) cos9e 11 ^. Variations with respect to A, p and r lead to the following set of 
equations of motion: 

= dn 2 p(uj) -1 , (2.23a) 

= ±+r 2 (u)-2Xtf(z) 

- -^ /d0 3 fdn' 5 p(Lj')log [f 2 (oj) + r 2 (uj'f - 2 f(w)f(a/) cos <p] , (2.23b) 

= Hu) _ 4# - A [™ 3 /dm f(^M-y)^^] . (2 . 23c) 

v ; 0r(w) 47r 4 7 6 J b f 2 {oj) + f 2 {oj')- 2r{uj)f{uj') cos if y ' 

It is far from obvious how to deal with this system of integral equations in the general case. Take 
notice that the equations are non-linear in r, and entangle p and r in a non-trivial fashion. 

6 



Nevertheless, for the ground state wavefunction where / = 0, it is known that the functional 
fflf is extremised by |7| 

A = iV 2 A = Ar2(io g Ar + ^_iog2) , p(u) = p = ^ , r(w) = r = ^ y . (2.24) 
This exact solution will play a key role in our subsequent discussion. 

3. Analytical solutions 

In this section, we will present analytical solutions to (2.23) for / ^ and hence, the most likely 
distributions of eigenvalues for excited wavefunctions ipf (i.e. eigenfunctions of the Hamiltonian 
H). We shall also compare our analytical results against numerical calculations. 

3.1. Perturbative expansions 

The starting point for our considerations is the ground state configuration (2.24): We consider 
wavefunctions that can be regarded as a slight perturbation of the ground state. For them we 
engineer a solution to the equations of motion as a perturbative expansion around the ground 
state solution. For this expansion to be consistent, we must assume that the function / = f(z) is 
'small' when evaluated on solutions to (2.23), i.e. we introduce some small parameter e and write 

/(*) := er 2 F{z) => \f(z)\ « r 2 = £ , (3.1) 

with F being holomorphic in z. The unknowns A, p and r are then expanded in powers of e 
according to: 

A = N 2 (A + eA 1 +e 2 A 2 + ---) , 

p(u) = PQ {l + e Pl {uj) + e 2 p2{uj) + ---) , (3.2) 

r(w) = r (l + eri(w) + e 2 r 2 (uj) -\ ), 

where Ao, po and r$ were given in (2.24). 

The equations of motion (2.23) can now be solved order by order in powers of e to eventually 
arrive at a perturbative solution of the form (3.2). The upshot of this expansion is that at any 
given order k + 1 (with k > 0), the integral equations for pt+i = pk+xioS) and r^+i = rk+i(w) 
are decoupled and linearised. Specifically, upon substituting (3.2) into (2.23), a few algebraic 
manipulations show that 

= dn 2 Pk+i(u) , (3.3a) 

Pk+i(to) = Afc+i - dn' 2 Ki(u},u}')p k+1 (uj') , (3.3b) 

Rk+i(co) = -r k+ i(uj) + dn' 2 K u (uj,uj')r k+1 (u}') . (3.3c) 
Here, we have introduced the integral kernels 

Ki(u,u') := — -5 /d^ 3 /d^3 log(l -cosy) , 

K u (oj,oj') := -L [anJdn'z—*: . 

27T° 7 7 1 — cos y 



(3.4) 



The remainders Pk+i = Pk+ii^) and Rk+i = Rk+i(w) are functions which depend only on the 
solutions A/, pi and 77, with I < k. They can be obtained order by order and the first few 
expressions are 

Pi(u) = MtF(z) 



1 



#iH 



(If 



+ - I dtt 2 Ri(oj) , 

<ReF(z) 



(3.5) 



r=r 



and 



ft(w) 



ife(w) 



J r\{u) + - / dft 2 n(w) [2pi(w) - ri(w)] 

O 7T 



+ 



^ fdn' 2 K ll (u;,io')r 2 1 (u') + ^ I dtt 2 R 2 (uj) , 



n(w) 



<9f 



5KeF(j8) - -ri(w) 2 + ri(w) / dO' 2 tf n (w,u/) pi (a/ 



(3.6) 



1 



dQ 2 K n(^,^')ri(^') [2pi(w') -n(u/)] . 



Notice that upon integrating equation (3.3b), the coefficients Afc + i for the Lagrange multiplier A 
are given by 



A 



fc+i 



7T 



d£l 2 P k+1 (uj) , 



(3.7) 



since JdQ 2 Kj(uj, a/) = JdQ 2 Kj(to, a/) = const and JdQ< 2 pk+\{oj) = 0; see Appendix A for more 
details. Hence, at any given order k + 1, the coefficient Afe+i is determined by the solutions A/, 
pi and 77, with I < k and does not depend on pu+i and rk+i- 

Equations (3.3b) and (3.3c) are Fredholm integral equations of the first and second kind, 
respectively. Such kind of integral equations can in principle be solved if the eigenfunctions of 
their kernels are known. We derive the eigenfunctions of Ki(u,u/) and Ku(uj,uj') in Appendix 
|A~1 where we also collect some facts about these kernels. Here we just give the results of these 
considerations. It happens to be that 



for m E Z and a E No with 

9 m , o (0) := cosl m ^Pi 1 'l m D(cos20) and $ m (0) := exp( 



lmc 



(3.8a) 



(3.8b) 



where Pa are the Jacobi polynomials (see Appendix B for definitions and properties), are the 
eigenfunctions of both kernels, though with different eigenvalues: 



dQ 2 Ki t n(uj,(J) i B m>a (u)') = \™ii^ m ,a(uj) , 



(3.9) 



with 



12 
-24 



log 2 
(\m\ +2a-l)\ 



Al' 



m\ +2a + 4)! 
m\ +2a)\ 



for (m,a) = (0,0) 
for (m,a) / (0,0) 



(3.10) 



11 (\m\ +2a + 3)! ' 

Notice that the functions (3.8) form a complete orthogonal basis for functions defined on the 



hemisphere given by uj = (6, (/>); see Figure 1. 



Therefore, we may expand the remainders appearing in (3.3) in terms of (3.8), 

Pk+i(w) = y^p^i^m,a{ u ) and %iM = y^Ci^M ■ ( 3 - n ) 

Since Pk+i is real, we must have Pj™?i = (P^j!^'")* and similarly for Rk+i] '*' indicates complex 
conjugation. Likewise, A& + i, p^+i and r^+i can be expanded in terms of (3.8). Upon replacing 
all of these expansions in the integral equations (3.3), we find that Afe+i, Pk+i and r^+i are given 
by 

Afc+i = P° k t , (3.12a) 



-.rii.a 



P 

Pk+M = ~ Y, ^ $ «»». (3.12b) 

(m,a)^(0,0) ! 

rfc+iCw) = ^y-^^H. (3.12c) 

m,a 3 + A II 

It should be stressed that in general there is no guarantee that the expressions (3.12b) and (3.12c) 
would define square-integrable functions. For that to happen, the functions Pfc+i and Rk+i should 
meet certain criteria. In particular, (3.3b) is a Fredholm integral equation of the first kind with 
symmetric kernel and a complete set of eigenfunctions. Therefore, certain existence theorems 
apply. Specifically, (3.3b) has a unique /^-solution given by (3.12b) if and only if the infinite 
series 



pm,a ^ 
"k+1 



\m,a 



(3.13) 



(m,a)^(0,0) 

is convergent. For more details, see e.g. |15|. 

In summary, provided certain criteria are met, we are able to solve the saddle point equations 
(2.23) analytically by means of the perturbative expansions (3.2). 

3.2. Example A: Monomial deformations 

Let us now consider an example and take f(z) oc z p for peN. Firstly, these are the simplest 
distortions of the ground state wavefunction. Secondly, and more interestingly, these wavefunc- 
tions should correspond to simple gravitational duals. Indeed, they are expected to be in cor- 
respondence with the LLM geometries obtained from a simply connected droplet possessing a 
non- vanishing p harmonic moment |10|. 

For this to constitute a small perturbation, f(z) should be of order erg, with £<1. Thus, 
the constant of the monomial is taken to scale with N. In particular, we consider 

f(z) = erl(— J => F(z) = (—} and <ReF(z) = (— J cos p 9 cos pcj) . (3.14) 

Notice that since 9KtF(z) is an even function in <f> and since Kjji(8, 4>, 8' ', (/)') = Kjji(0, —(f), 8', —</>') 
as a direct consequence of (2.21) and (3.4), the functions p and r will be even functions in 4>. 
Hence, the series expansions (3.12) reduce to expansions in terms of ^\t^ m)a with m > and real 
coefficients (i.e. we only need the cosm</> terms in the Fourier expansion). 
Therefore, the expressions (3.5) for the remainders Pi and R\ are given by 

Pi(w) = cos p 8cospcj) = ^te^oM , 

(3.15) 
Ri(ui) = pcos p 8cosp<p = p9\ety p .o(oj) . 



Using these results and the series expansions (3.12), we then obtain 

Ai = 0, 

1 (p + 4)! 



Pi(w) 

n(w) 



24(p-l)l 

3p(p + 3)! 
2(p + 3)! + 24p! 



cos r t/ COS 1 



cos r t/ cos ; 



(3.16a) 
(3.16b) 
(3.16c) 



2p,0 
P /»2 



0,0 0,1 

P 2 P2 



0,2 
92 



0,3 
/>2 



0,4 
P2 



2 
3 



925 
3 

435375 



121 

a 29458800 
^ 1369 









10 



1S5 



2430 89910 
121 





87075 q 



847 847 
40320 7560 481600 233800 



1369 37 



1369 



1369 



Table 1: Coefficients of p 2 (co). 



p 


2p,0 

r/' 


0,0 
r 2 ' 


0,1 

r 2 ' 


0,2 

r 2 ' 


0,3 

r 2 


0,4 
r 2 ' 


2 


2875 
592 


115 
32 





575 

1184 








3 


42795 
1892 


12681 
968 


1125 

484 


492075 
250712 


8559 
13244 





4 


14391825 
228623 


136885 
4107 


43840 
4107 


348705 
101306 


477680 
176601 


685325 
1371738 



Table 2: Coefficients of r 2 {oS). 

Next we would like to compute the solution to order e 2 . As it is clear from (3.6) and (3.16), 
in order to compute A2, p 2 and r 2 we need to expand (cos p 6 cos p4>) 2 = [fRt^f p fl(uj)] 2 in terms of 
the eigenfunctions (3.8). We find 



9 =0 



(p-q)Kp + q + 2)\ 



Therefore, using (3.6), the non-vanishing coefficients p™' a and R 2 ' a , 

P 2 {lo) = YsK^^mA^) and R 2 (oj) = ^i?™' a <He* m , a (G 



m,a 



are computed to be 
bpCpfl 



hi (2 



V + 



Att a 



11 I "p 






2p,0 



p°P,9 



+ A 



0,g 



i?o 



OpCp,q 



pfl 



— - b p + OpAjj 
P 2 - o ^p + a P Ag c 



for (m,o) = (0,0) , 

for (m, a) = (2p, 0) , 

for (m, a) = (0, 1 < q < p) , 

for (m, a) = (2p, 0) , 

for (m, a) = (0, < q < p) 



>V \ \2p,0 

2"J Al1 



p 



Mi 



(3.17) 



(3.18) 



(3.19a) 
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where 



1 (p + 4)l = 3p(p + 3)! _ (g + l) 2 p! 2 

24(p-l)!' p '"" 2(p + 3)\ + 24p\ &U Cp ' q ''~ (p - q)l(p + q + 2) 



and C VA ■= ,„ Z,J , , oM ( 3 ' 19b ) 



and the eigenvalues A^'" were given in (3.10). From these expressions, we obtain 

A 2 = p°.°, p™>° = _ 2 for (m;fl) _, (0>0) and r m,a = _^_ ^ 

as the only non-vanishing coefficients appearing in 

p 2 (co) = Y,P2 ,a ^ m ,a(u) and r 2 (w) = ^r™' a $He* m , a (a;) . (3.21) 

m,a m,a 

We refrain from writing down these expressions explicitly, as they are lengthy and not particularly 
illuminating. Instead, we only list the particular values for p = 2,3,4 in Tables 1 and 2, 

In what follows, in order to compare against the analytical results, we will compute these 
coefficients out of numerical solutions. To have a numerical handle on the problem, we can 
retrace our steps from the continuum limit to finite N. For finite N, the extrema of Jiff will 
not dominate the partition function (2.11a) anymore as the saddle point approximation is only 
valid in the thermodynamic limit iV — > oo. Instead, for a finite TV" description, one should take 
into account all possible configurations weighted by their respective probability. In order to do 
that in a good approximation, Monte Carlo simulations with a Metropolis criterion can be used 
to simulate the distribution of probabilities: Starting off from a random distribution of particles, 
a new one is generated by a small random perturbation. The new configuration is accepted if 
e~ f (where Jiff is the discrete Hamiltonian (2.11b)) is larger than a random number in the 
interval [0, 1] and rejected otherwise. By iterating this algorithm a large number of times, a 
typical (i.e. most likely) configuration of particles will be obtained at the end. This method was 
used to describe different wavefunctions of the matrix model we are dealing with for finite N in 
works by Berenstein and collaborators |6,12,13|. 

At this point, we should emphasise that the remainder of this section is not aimed at finding the 
finite A" description of wavefunctions. Instead, we will use the finite A^ problem as a discretisation 
of the equations of motions of the continuum limit, which, of course, would be valid in the A" — > oo 
limit. In other words, we will consider the finite A" Hamiltonians Jiff and still look for their 
extrema numerically and regard these as numerical approximations of the perturbative analytical 
solutions presented previously. To do that, we will simply modify the Metropolis criterion to 
accept new configurations only if 5 Jiff < 0. In this way, the iteration procedure will produce 
configurations approaching the extrema of Jiff. All the simulations in this work that use this 
modified Metropolis criterion were performed in Mathematica. 

Using that method we obtained numerical approximations for the cases p = 2,3,4 with e = 
0.01. In these three cases, we discretised the densities using A" = 1000. For such small e, the 
main dependence of f will be that of the linear order in e. If one plots f versus cos p 9 cos p4>, 
the points should approximately lie in a line. In Figure 2 we produce those plots on top of their 
linear fittings. The agreement of these numerical fittings with the analytical coefficients is very 
good. Certainly, beyond the the leading approximation, r is not a function of cos p 9 cos p4> only. 
By fitting f as an expansion of the eigenfunctions kicking in at higher orders in e one can obtain 
the corresponding numerical coefficients. For instance, for p = 2, we get the coefficients 

(r$>\r°/,r /,r°f)\_ lcal « (5.1,-3.7,0.02,0.58) (3.22a) 
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cos 2 9 cos 2<f> 



cos J 9 cos 3<j> 



cos y cos ' 



Case p = 2: 

This is very close to the analytical solu- 
tion v/500(l - ^e 2 - §e cos 2 (9 cos 20) « 
22.35 -0.559 cos 2 cos 20. 



Case p = 3: 

This is very close to the analytical solu- 
tion y/500(l-^±e 2 -§e cos 3 9 cos 30) « 
22.33 -0.915 cos 3 cos 30. 



Case p = 4: 

This is still close to the analytical solution 
V / 500(l " "^ " W £ cos 4 cos 40) « 
22.29 - 1.269 cos 4 cos 40. The slight de- 
viation is due to the fact that for larger 
p the numerical value of the (neglected) 
higher-order coefficients is larger. 



Figure 2: Numerical solutions for p = 2,3,4 and e = 0.01 obtained with N = 1000. The solid 
curves represent linear fittings of the data points. 



which are fairly close to the coefficients in first row of Table 2 : 

(rr,rr,rr,r°/)L, w = (»"W.0,ffi - (4.7,-3.6,0,0.49) . (3.22b) 

It is also possible to compute the density coefficients for the numerical solutions using 5 

JV 



P 



m,a 



1 



with the normalisation constant 



dQ 2 /5(w)^e^ m , a (w 



1 1 



' i=l 

ir (a + l)(l + <5 m0 ) 



5^JHe* m>a (wi) , (3.23a) 



2 (2a + m + 2)(a + m+l) ' 
5 The p m ' a s are the coefficients of the full p, i.e. p{u) = E m a p m ' a 9^* m ,„H 



(3.23b) 
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where 5 is the Kronecker delta. For p = 2, we obtain 

(pl' ,pt'° ,P2'\P2' 2 )L mcrical « (30.1,306,10.1,30.0) (3.24a) 

in good agreement with the coefficients in Table 1: 

(p?'°,pf,P^,P^)L nBlytical = (30,¥,10,^)« (30,308.3,10,30.8). (3.24b) 

3.3. Example B: Logarithmic deformations 

Let us now consider a wavefunction with a logarithmic potential f(z) = Qlogz. In this case, 
the wavefunction should correspond to the BPS operator det-Z^ |2|. In turn, this operator is 
believed to be dual to the annular LLM geometry whose droplet inner and outer radii are \/Q 
and \/N + Q, respectively. 

One could expect this case to be simpler because the (^-dependence drops out of the problem 
(notice that the real part of / is independent of <p) and the functions f and p will only depend on 
0. On the other hand, the distribution cannot be regarded as a small distortion of the spherical 
for all values of the angle (since / has singularities) and therefore we cannot use (3.2) to solve 
the equations perturbatively. 

Nevertheless, there is a simplification since for a logarithmic potential, the term —n^- appear- 
ing in the 'radial' integral equation (2.23c) does not depend on 9 explicitly. Therefore, (2.23c) 
can be solved exactly by assuming a constant radius r: 



N 
r = r Q = Yy+Q- (3.25) 

To arrive at this result, we have used the constraint (2.23a) to eliminate p. We should emphasise 
that this solution is an exact solution for any value of Q. 

Pleasingly, this appears to be consistent with Monte Carlo simulations for finite values of 
N |6,13|. In Figure 3, we have depicted R$q := \Jx\ + x\ versus -R1234 := \Jx\ + x\ + x\ + x\ 
for N = 2000 particles and Q = 20. The plot on the left corresponds to a typical (i.e. most 
likely) configuration obtained from a Metropolis algorithm (as done in [6, 13|). The solid curve 
corresponds to -Rf 2 34 + -^56 = r Q w ith Tq = 1020. Of course, the agreement would improve 
for simulations with larger N. As before, one can modify the Metropolis criterion to accept new 
configurations only if 5J#f < 0. In this way, the iteration will produce configurations approaching 
the extrema of Jtiff . As shown in right plot in Figure 3 , these configurations fit very well the large- 
A^ value of the constant radius (3.25), 6 

Now that we have solved exactly the radial integral equation with (3.25), we are left with a 
linear integral equation for the density p(9), 

A- ^log(l + cos2fl)-- / d0'cos0'sin 3 / .K'j(0,0 / )/5(0 / ) = 0, (3.26a) 

A 2 Jo 

with 

A := A2+^-f lo gy- lo g 2r Q and K I (9,9') := f fyKjfau/). (3.26b) 

J 
6 In fact, several simulations for different TV and Q values show that the particles of extremal configurations sit 

at tq = yJH=±+Q. 
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^1234 




*" fil234 



(b) 



Figure 3: On the left, a typical configuration for N = 2000 and Q = 20 obtained by using 
a Metropolis criterion in the numerical simulation. On the right, a configuration extremising 
Jiff which is also obtained for iV = 2000 and Q = 20. In both pictures, the solid curve is 
1020. 



f? 2 

-"■1234 



+ ^56 



Using the expansion (A. 15) of the kernel Kj, we may re- write the kernel Kj as an infinite series 

oo 

K!{0,ff) = 4^A?' a (l + a)Pi 1 '°)(cos2^)P ( ( 1 ' )(cos2^) , (3.27) 



a=0 



>0,a 



where A^' was given in (3.10). Likewise, we could integrate (A. 6) to arrive at an infinite series 
in terms of hypergeometric functions. 

Equation (3.26a) is again a Fredholm integral of the first kind. In general, the range of validity 
of such an equation does not need to coincide with the interval of integration. In fact, (3.26a) 
cannot hold for all < 9 < § as can be seen as follows: 7 Firstly, (3.27) implies that Kj is 
bounded 8 since with |P„ ' (cos 29) \ < 1 + a for < 9 < § we find 



\Kj(e,tf)\ < 4X>?' a |(l + a) 3 = 41og2 



(3.28) 



a=0 



Secondly, the integral (recall that p > 0) 



tt/2 



d9' cos 0' sin 3 6' K T (e, 6') p(0') 



< 



tt/2 



&6' cos 0' sin 3 & | Ki{B,&) \p{&) 



< 41oe2 



tt/2 



d^cos^sin 3 e / / 5(^ 



7T 



log 2 (3.29) 



is also bounded for all < 9 < \. In the last step of this derivation, we have used the normalisation 
of the density p. However, the logarithmic term log(l+cos 26) appearing in (3.26a) is not bounded. 
Therefore, we conclude that (3.26a) cannot be solved for all < 9 < ^. Note, however, that this 



7 To arrive at this conclusion, one may also argue differently. Upon expanding log(l + cos 20) and p in terms 
of the Jacobi polynomials Pi ' (cos 20) and using (3.27), one quickly realises that the obtained series p{9) — 
^2 paPa ' (cos 29) is not convergent as the coefficients p a do not satisfy the criterion (3.13). Hence, no solution 
exists satisfying the equation for all < 9 < ^. 

8 Alternatively, one may deduce this directly from the expressions (A.6). 
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does not contradict the intuition that a function p extremising JfCf should exist. Rather, since 
p > 0, the extrema of J#f could lie in the boundary of the space of the allowed configurations 
p. Indeed, the above analysis indicates that the extrema of J#f must be found in the boundary 
of the configuration space, i.e. the density p has to vanish in some region U C [0, |] and in this 
region U, the integral equation will not hold. 

This is consistent with the numerical simulations, where p appears to be vanishing in the region 
Ue = {9 | 6q < 9 < ?} for some 9q. Altogether, we are left with the following (one-dimensional) 
integral equation for p: 

o ir r e ° 

A-^log(l + cos2fl)-- / d9' cos 9' sin 3 # Ki(0, 9') p{9') = 0, for < 9 < 9 (3.30) 
A 2 J 

for some value 9q = #o(]v) which one has to determine consistently together with p. 9 To find 9q 
and p, one should solve (3.30) for generic 9q ( n °t for all 9q a solution should be found). Among 
those solutions, one then should look for the one that extremises Jtff. 

Unfortunately, solving the Fredholm integral equation (3.30) for arbitrary 9q is a very difficult 
problem. So far, we have not been able solve the eigenvalue problem of the kernel (3.26b) for any 
other value than 9q = % . For that reason, we shall determine p numerically in the following. 

Let us first make a change of coordinates according to x := cos 29, so that the normalisation 
and the integral equation (3.30) read as (xo := cos26*o): 

1 4 
dx (1 — x) p(x) , (3.31a) 

= A - % log(l + x) - ^ / dx' (1 - x') K I (x, x') p(x') . (3.31b) 

N 16 7 xo 

These equations can be obtained from variations of 10 

jf t = a (J\ a - x) m - 1) - 1 [mi - ,) io g( i + ,)«,) 

dx (1 — x) I dx (1 — x ) Kj(x, x) p(x) p(x') . 



32 



(3.32) 



Note that the full range — 1 < x < 1 is used in Jtff. The restriction to xq < x < 1 will be produced 
for configurations of p that exactly vanish for x < Xq. 

To extremise Jif, we can discretise the problem by thinking of a large number L of 'particles' 
in the interval [—1,1] and trading back integrals into sums 



/_ 



i 4 L 

dx(l-x)p(x) -»• —Y, (3-33) 



and hence, 



&t = -^E l0 S( 1 + ^)-^l2 Y^K^xj). (3.34) 

To obtain the extrema of J^ ; we now perform Monte Carlo simulations with a modified Met- 
ropolis criterion (i.e. we only accept new configurations if SM'f < 0). The big advantage of 



9 Notice that #o(]y = 0) = f . Furthermore, a necessary condition is that A — S log(l + cos2^o)| < log2 as 
follows from (3.29). 
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Notice that Jiff is basically a redefinition of Jiff . 
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(l-x)p(x) 



(l-x)p(x) 




1 



N 500tt 

A m -0.1103 
x & -0.8495 



0.0006 




(b) 



(l-x)p(x) 




0.0032 



(l-x)p(x) 
10 



4 J 






N 50ir 

Ass -0.1139 
x « -0.5803 



0.0064 



^/rfl 1 1 ITh I n rfTnTn rTrrfTrrn*~***r+ . _^_ _ .. j ^ 



(c) 



(d) 



0.0095 




fcWWVW /VS/\/W\/l I 



(l-x)p(x) 

I: 



Q 



1 



0.0127 



JV 25tt 

A« -0.1173 
x « -0.4566 

,— fTffTTI nHTfTVl-Tli>T--n-r-r>- n j , ^ 



0.5 



(f) 



Figure 4: Numerical solutions for (1 — x)p(x) for various values of ^. The number L of 'particles' 
for doing the numerics is always 2000. The solid line represents an interpolating function obtained 
by numerically expanding the solutions p in terms of the Jacobi polynomials. Notice that the 
analytical result for § = is p{x) = £ « 0.6366, A = ^ - log 2 sa -0.1098 and x = -1. The 
small oscillations are an artifact of the discretisation: They become smaller as L is taken larger. 
These results agree with those found previously in |6,13|. 

the numerical approach considered here when compared to the one for the original Hamiltonian 
(2.11b) is that we have effectively reduced the six-dimensional problem to a one-dimensional 
problem. However, there is also a slight disadvantage in that the interaction term appearing in 
Jtiff is more complicated. Notice that the numerical value for the Lagrange multiplier A (and thus 
A via (3.26b)) can be found from the discretised version of (3.31): 



= A-^log(l + a* 



1 L 



(3.35) 



j=i 



Notice also that this equation can be used as a criterion to estimate how good our numerical 
solutions are (since A should be constant for all Xi). In Figure 4, we have depicted our numerical 
results for various values of ^. The number L of 'particles' was chosen to be 2000. 
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4. Conclusions and outlook 

We have considered a particular quantum mechanics model of commuting matrices which is be- 
lieved to describe 1/8 BPS states in SU(iV) M = 4 SYM theory at strong coupling. The wavefunc- 
tions of this model are expected to be dual to type IIB supergravity solutions that asymptotically 
approach AdSs x S 5 . The probability densities for the wavefunctions are interpreted as partition 
functions of certain matrix models. Different wavefunctions correspond to different potentials in 
the associated matrix model Hamiltonian. 

Specifically, we focused on the large- N or thermodynamic limit in which the probability dens- 
ities are dominated by the saddle points of the matrix model. Then, we have solved analytically 
the saddle point equations for a family of wavefunctions that are in correspondence with specific 
LLM geometries. 

The starting point of our consideration was the ansatz (2.17), where it is assumed that the 
probability densities are supported on particular hypersurfaces in R 6 . We then constructed a 
perturbative approach that allows for an analytical treatment of the saddle point equations. We 
used monomial potentials to illustrate how the perturbative method works. By perturbing around 
the ground state solution, we constructed the solution up to and including second order in the 
deformation parameter. 

We then considered a logarithmic potential, which does not admit a perturbative solution 
to its saddle point equations. In spite of that, we could also obtain partial analytical results 
in this case. We found an exact solution (3.25), which holds true for any value of Q, for the 
radial function of the ansatz (2.17) and we reduced the problem to a linear integral equation for 
the density p (3.30). To provide a full analytical solution in the logarithmic case, one should 
diagonalise the integral kernel given in (3.26b) for a generic interval < 9 < 9q with 9q < f . 
Finally, we have compared all our analytical results against numerical simulations and found very 
good agreement. 

The main issue which we have not discussed here but we hope to report on in the future 
concerns the extraction of geometry from our analytic solutions. Recall that for the ground state 
solution, the saddle point configuration is given by a uniform probability distribution supported 
on a five-sphere in 1R 6 and an explicit geometrical meaning to this five-sphere was given in |2,7|. 
The generalization of this identification is not that obvious for the excited wavefunctions. To begin 
with, the dual LLM geometries do not have a factorised five-dimensional compact geometry. Put 
it differently, the compact factor is different for different sections of the LLM geometry. At 
least one would expect to be able to reproduce the geometry of the section corresponding to the 
LLM plane. Then, there is also the question of how to read geometry out of the distribution of 
eigenvalues. There is obviously the hypersurface where the density is supported, but the density 
function must also enter somehow. For instance, although the supporting hypersurface in the 
logarithmic is a five-sphere, the compact factor in the LLM plane of the annular LLM solution 
ought to be something else. 

It would also be interesting to extend our perturbative method to solve analytically the saddle 
point equations of more general wavefunctions. For instance, potentials that depend also on other 
holomorphic coordinates are supposed to describe wavefunctions of 1/4 and 1/8 BPS states. 
This would shed light on the AdS/CFT dictionary for those cases (see e.g. |16| and references 
therein). In addition, it would be very interesting to use our perturbative method to characterise 
excited wavefunctions in other matrix quantum mechanical systems |4,5|, which according to the 
AdS/CFT correspondence are dual to supergravity solutions with other asymptotics (e.g. AdSs x 
T 1 ' 1 ). 
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Appendices 



A Eigenvalues and eigenfunctions 

In this appendix, we shall derive the eigenvalues and eigenfunctions of the kernels K\ and K\\ 
introduced in (3.4). Before tackling the problem, we begin with some preliminary considerations. 

Preliminaries 

In the definition (3.4) of the kernels K\ and Kn, we are integrating over (a, £i, £2) and (a' , £[, £ 2 ) 5 
i.e. over all possible positions of two points on a three-sphere. The dependence of (p on them is 
through the relative angle between two points on the three-sphere, so we can fix one of the points 
arbitrarily (say a' = and £[ = 0) and multiply by the volume, 2ir 2 , of the three-sphere. In the 
mentioned choice, also the integral over £2 becomes trivial and we obtain 

o / i7r /2 / 1 2-7T 

Ki(oj,oj') = — 2 / da d£i sinacosa log(l — cosi^o) , 

71 Jo , Jo ' (A.la) 

9 rn/2 ,-2-k i v ' 

Ku(u),u)') = — k / da d^isinacosa 

7T Z J Jo 1 - COS ip 

where 

cost/?o := cos cos 0' cos (0 — (j)') + sin # sin #' cos a cos £1 . (A. lb) 

There are two (equivalent) ways of computing (A.la). Firstly, we may expand the expressions 
log(l — cos(/?o) and 1/(1 — cos^o) in powers of cos</?o: 



00 1 00 

Ki(uj,uj') = — > — Kj(co,u)') and Kji(uj,U)') = > Kj(uj,lo') , (A. 2a) 

with 



j=\ J j=0 



fKJ2 p2ir 



A short calculation reveals that 



Kj(cj,uj') := — s- / da d£i sin a cos a cos- 5 ' tpo ■ (A. 2b) 

^ Jo Jo 

zeals that 

tj- t i\ 2 . . .... • (\ — j —j tan 2 9 tan 2 9' \ . . 

^■(w.w') = - cos 6> cos #' cos </>-<// P 2^1 ^HS^f; 2; , (A.3 

7T \ 2 2 cos z (<p — <p ) J 

where 2^1 is the hypergeometric function. Notice that either the first or the second argument 
of the hypergeometric function appearing in (A.3) is a non-positive integer. This implies that 
2.F1 (-9-^) -%-'■> 2; •) is actually a polynomial in the last argument. 

As an alternative to the above expansions, one may perform the integrals in (A.la) directly. 
In particular, they can be re- written as 

o ^1 p2n 

Ki(uj,uj') = -2 / dt / d£itlog(A-BtcosZi) , 
Xii(w,u/) 



^ 2 . 


rl /*27T 
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- Bt cos £1 



where 



A = A(uj,uj') := 1 — cosOcosO' cos(cf) — (f)') and B = B(uj,uj') := sin 9 sin 9' . 
Notice that \B/A\ < 1 for all (a;, a;'). Using (|a| < 1) 



/ dx log(l — acosx) = 27rlog 2(1 + v 1 — a 2 ) , 
Jo 



2n 

dx 



1 



2tt 



/ 1- acosx y/i - a 2 

the kernels (A.4a) are given by 

1 -(a + VA 2 - B 2 ) - I + 4o ( A - VA 2 - B 2 



Ki(u,u') 



2 



log- 



2 B 2 



, i S 2 / 3 £ 2 \ 



7T 



where 3.F2 is a generalised hypergeometric function. 

Notice that as a by-product we have obtained the non-trivial identities 



•E^^lf^,^^:// 2 ) - lo P ,(.l-,-) 






2 ' 2 

1 - i - 3 



x v I 3 xv 



(A.4b) 



;2;y' 



2 ' 2 
between hypergeometric functions. 



1 _ (1 , xV 

h 2-fl o' 1 ' 2 ' 7l ^2 

I — x V 2 1 — xY 



(A.5) 



(A.6) 



(A.7) 



Eigenvalue problem 

Let us now derive the eigenvalues and eigenfunctions for the kernels K\ and K\\. To this end, the 
form (A.2) turns out to be very adequate, since the sub-kernels (A.3) can actually be diagonalised 
simultaneously 11 In fact, 



4Q! 2 Kj{u,(J)* m <a{(J) 



(1 + (_1)M+J) 2 -Jr(j + 1) 

r(|-M_ a + 1 )r(| + M + a + 3 ) 



* m ,aM , (A. 



=: A" 



where 



* m , a ( w ) = cos |m| 0Pi 1,|m|) (cos20)exp( 



lmc 



(A.8b) 



>(«,/») 



for to € Z and a G No- Here, Pa ' are the Jacobi polynomials (see Appendix B) and V denotes 
the Gamma function. 

"Therefore, any kernel of the form ^j. cijKj with some constant coefficients a,- can be diagonalised in this 
manner. 
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These eigenfunctions form a complete orthogonal basis for functions defined on the hemisphere 
given by u; see Figure 1. Indeed, upon introducing the delta function 

S(u,u/) := l^] -J{9 -9')8{(t>- if/), with - [an' 2 5(oj,oj') = 1, (A.9) 

2 sin v cos v vr J 

the completeness relation is given by 

E(\m\ + 2a + 2)(|m| + a + 1) Al , . 

^— ! 2a + 2 -*m,aM[*m,a(w')] = %,"')> (A.10) 

where '*' denotes complex conjugation. This follows from the completeness relation (B.21) of the 
Jacobi polynomials. Furthermore, the orthogonality relation reads as 

I J™ 2 * m , a M[*n, b M]* = (H+2a ^, 2 m|+Q+1) W- , (A.11) 

where 5 is the Kronecker delta. This expression can be obtained from the orthogonality relation 
(B.20) of the Jacobi polynomials. Altogether, we may obtain an orthonormal basis by setting 

. n , /v |m| +2a + 2)(|m| +o + l) 



tt m , a (a;) H- *^ a (w) := A /^-! ^LJ ^* m , a (a;) (A.12) 



for which we have 



E f ^( w )[ f m,<.( w ')]* = ^(^^') and ljdn 2 ^ a (u,)[V ntb ( U )]* = 6 mn 8 ab . (A.13) 



Next we would like to derive the eigenvalues Ajj". They are given by summing up the 
m 
•J 



eigenvalues A" 1 '" appearing in (A. 8): 



oo 

V - " 2^ 7 A i - 1 



,=1 J 



— -log 2 for (m,a) = (0,0) , 

(|m|+2a-l)! for 
(|m|+2a + 4)! v ; r v ; (A.14) 



OO 






(|m| +2a)! 



(|m|+2a + 3)! ' 

Using the eigenvalues A] 7 ^'" and the normalised eigenfunctions ^^a, the integral kernels K\ and 
K\\ may be expressed as 

Kmfau/) = ^£A™^>) [<,>')]*• (A-15) 

m,a 

Notice that upon using (A.6), one easily computes the traces tr K\ and tr K\\. They can also 
be obtained by summing up the eigenvalues (A.14): 

trAT.,11 = fdShKjpfau) = ^Ayf = | _1 ° g2 ' (A.16) 
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B Jacobi polynomials 

The Jacobi polynomials Pa for a, (3 € R (with a, (3 > —1) and a e No are solutions to the 
ordinary differential equation 

(1 - x 2 )y" + ((3-a-(a + f3 + 2)x)y + a(a + a + f3 + l)y = (B.17) 

and they can be obtained from the hypergeometric function 2-F1 according to 

pK«(x) = (a ^, 1)a 2F1 (-a, l + a + P + a;a + l; ^-\ , (B.18) 

where (a) n := T(a + n)/T(a) is the Pochhammer symbol. We are particularly interested in the 
polynomials Pa which are: 

P^ m \x) = 1 , 



P^ m \x) = 2 + i(3 + m)(x-l), 

P^' m \x) = 3 + |(4 + m)(x-l) + |(4 + m)(5 + m)(x-l) 2 
P 3 (1 ' m) (x) = 4 + 3(5 + m)(x-l) + i(5 + m)(6 + m)(x-l) 2 
+ ffi(5 + m)(6 + m)(7 + m)(x - l) 3 , 

Furthermore, the Jacobi polynomials form an orthogonal basis with 

2 a+ ? +1 (a + 1) Q 



(B.19) 



/1 
dx(l-x) a (l + ^P a (a ^ ) (x)P b (Q ' /9) (x) 



<U . (B.20) 



2a + a + /3 + l(/3 + a + l) 
The completeness relation is given by 

- {a+p+2 a +1){P+ a + ru p ,„„ <„„ = 2"*" J(l _ , (B . 21) 

^ (tt+l) a a \ ) a \yj (l_ rE )a(l +x )/3 V V) \ I 

for |x| < 1 and \y\ < 1. 

Finally, we record the following useful relations: 

pM(- x ) = (-l) a P^(x) (B.22) 

and 
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